Fluctuation dynamo based on magnetic reconnections 
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We develop a new model of the fluctuation dynamo in which the magnetic field is confined to thin flux ropes 
advected by a multi-scale flow which models turbulence. Magnetic dissipation occurs only via reconnections 
of flux ropes. The model is particularly suitable for rarefied plasma, such as the Solar corona or galactic halos. 
We investigate the kinetic energy release into heat, mediated by dynamo action, both in our model and by 
solving the induction equation with the same flow. We find that the flux rope dynamo is more than an order of 
magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic 
energy released during reconnections has a power-law form with the slope —3, consistent with the Solar corona 
heating by nanoflares. We also present a nonlinear extension of the model. This shows that a plausible saturation 
mechanism of the fluctuation dynamo is the suppression of turbulent magnetic diffusivity, due to suppression of 
random stretching at the location of the flux ropes. We confirm that the probability distribution function of the 
magnetic line curvature has a power-law form suggested by Schekochihin et al. 12511 . We argue, however, using 
our results that this does not imply a persistent folded structure of magnetic field, at least in the nonlinear stage. 



I. INTRODUCTION 

Dynamo action, i.e., the amplification of magnetic field by 
the motion of an electrically conducting fluid (plasma), is the 
most likely explanation for the omnipresence of astrophysical 
magnetic fields. Ohmic dissipation, however small, is essen- 
tial in order to achieve the development of the dynamo eigen- 
solutions and to smooth out the spatial variations of the mag- 
netic field. The evolution of the magnetic field B embedded 
in a velocity field u is governed by the following closed equa- 
tion: 

<9B 

— = Vx (uxB)+£B, (1) 

where C is an operator describing the magnetic dissipation. 

In rarefied astrophysical plasmas, such as the Solar coro- 
na, hot gas in spiral and elliptical galaxies, galactic and ac- 
cretion disc halos, and laboratory plasmas, an important (or 
even dominant) mechanism for the dissipation of magnetic 
field is the reconnection of magnetic lines rather than mag- 
netic diffusion l20ll . the latter modelled with C = 77V 2 (if 
i] = const). Discussions of astrophysical dynamos often refer 
to magnetic reconnection, but attempts to include features spe- 
cific of magnetic reconnection into dynamo models are very 
rare [see however Ht]. On the other hand, theories of mag- 
netic reconnection (and the resulting estimates of the plasma 
heating rate) rarely, if ever, refer to the dynamo action as the 
widespread mechanism maintaining magnetic fields. This pa- 
per attempts to bridge the gap between the two major top- 
ics of astrophysical magnetohydrodynamics (dynamos and re- 
connections) by developing a dynamo model which explicitly 
incorporates magnetic reconnections. 

The nature of the dissipation mechanism is important for 
the dynamo action as it affects the growth time of magnetic 
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FIG. 1: A schematic representation of the magnetic dissipation op- 
erator £ in Fourier space: usual diffusion Ck oc k 2 (dash-dotted), 
hyperdiffusion Ck oc fc 4 (dashes) and reconnections at a scale do as 
described by our model (solid). 



field, its spatial form and the rate of plasma heating by the 
electric currents. For example, dynamo action with hyperdif- 
fusion, L = —771V 4 (and with a helical u) has larger growth 
rate and stronger steady-state magnetic fields than a similar 
dynamo based on normal diffusion [5]. This is not surprising, 
as the hyperdiffusion operator, having the Fourier dependence 
of k A , rather than the k 2 dependence of normal diffusion, 
has weaker magnetic dissipation at larger scales as shown in 
FigH This allows the magnetic field to grow unimpeded by 
dissipation as magnetic dissipation is confined to relatively 
small regions. The release of magnetic energy in smaller re- 
gions (and larger current densities) in hyperdiffusive dynamos 
may also lead to a higher rate of conversion of kinetic energy 
to heat via magnetic energy. One of the aims of this paper is 
to demonstrate that this statement is especially true in the case 
of magnetic reconnections. 
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FIG. 2: (Colour online) The algorithm for inserting new trace parti- 
cles in a stretched (left to right) or contracting (right to left) magnetic 
flux tube. If the distance between any two trace particles (shown with 
red/open circles) exceeds a length scale d, a new particle is inserted 
between the two particles shown with a blue/filled circle. The label 
at each particle represents magnetic field strength at that location 
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FIG. 3: |B| at a specific position (y — 1.) in a shear flow 10 whose 
velocity is reversed at t ~ 20. The dotted line shows the analytic 
solution l|4j, and numerical results are shown with solid line. The 
initial field strength is B = 1. 



Magnetic hyperdiffusion also appears in the context of con- 
tinuous models of self-organised criticality in application to 
the heating of the Solar corona [6]. The aim of such models 
is to reproduce the observed frequency distribution of various 
flare energy diagnostics. As we show here, our model exhibits 
a power-law probability distribution of the magnetic energy 
release similar to that observed in the Solar corona. 

Magnetic reconnection may correspond to an even more 
extreme form of the dissipation operator than the hyperdif- 
fusion: here magnetic flux tubes dissipate their energy only 
when in close contact with each other, so that the Fourier 
transform of C should be negligible at all scales exceeding 
a certain reconnection length c?o (see Fig . [TJ . It is then natural 
to expect that dynamos based on reconnections (as opposed to 
those involving magnetic diffusion) will exhibit faster growth 
of magnetic field, more intermittent spatial distribution and 
stronger plasma heating. In this paper we consider dynamo 
action based on direct modelling of magnetic reconnections. 
For this purpose, we follow the evolution of individual closed 
magnetic loops in various flows (known to support dynamo 
action) and reconnect them directly whenever they come into 
a sufficiently close contact, with appropriate magnetic field 
directions. First results of our simulations can be found in |Q]]. 



II. THE FLUX ROPE MODEL 

We model the magnetic field by considering the evolution 
of thin flux tubes, frozen into a flow, each with constant mag- 
netic flux if). We first focus on the kinematic behaviour, where 
the velocity field is independent of magnetic field. Later we 
shall introduce the Lorentz force into the system to account for 
the back reaction of the magnetic field on the velocity field. To 
ensure that V • B = 0, we require that our flux tubes always 
take the form of closed loops. Numerically, we disctretize the 
loops into fluid particles and track their positions and relative 
order (i.e., magnetic field direction) by introducing a flag P, 
whose value increases along a given magnetic flux tube. Ini- 



tially the particles are set a small distance apart, 0.75c?, where 
d is a certain (small) constant length scale. If, during the evo- 
lution of the loops, the distance between two neighbouring 
fluid particles on a loop becomes larger than d, we introduce 
a new particle between them, as illustrated in Fig. [2] We use 
linear interpolation to place the new particle halfway between 
the old ones. For example if inserting a new particle x c be- 
tween particles x a and Xb the position of the new particle is 
given by, 

x c = i(x b +x a ). (2) 

The separation between the new particles is thus greater than 
0.5c? - a feature which will be important when we consider re- 
moving particles. The effective spatial resolution of our model 
is thus close to d. We shall discuss later a prescription for x c 
which is more accurate than Eq. @. 

Each particle is also assigned a flag B which denotes the 
strength of the magnetic field at that position on the loop. As- 
suming magnetic flux conservation and incompressibility, the 
magnetic field strength in the flux tube is proportional to its 
length. Initially the magnetic field is constant at all particles 
in each loop, B = 1. However, when a new particle is in- 
troduced, the magnetic field is doubled at certain particles, as 
shown in Fig. [2] Importantly, the field strength is increased at 
two out of three particles involved: this prescription emerged 
from our experimentation with various schemes, and allows 
us to reproduce the evolution of magnetic field strength in a 
shear flow. Conversely, when the flow reduces the separation 
of particles to less than 0.5c?, we remove a particle. The value 
of the magnetic field strength flag is also halved on the remain- 
ing particles in a manner consistent with the above algorithm. 
We have verified that this prescription reproduces accurately 
an exact solution of the induction equation for a simple shear 
flow. 

Results presented below have been obtained with a typical 
number of trace particles of order 10 4 . 
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FIG. 4: (Colour online) The shape of the flux tube stretched by the 
flow, given in Eq. l[3j at t — 2.7. Colour coding shows the magnetic 
field strength according the key (right), Bq = 1. 




FIG. 5: (Colour online) Reconnection occurs when the distance be- 
tween two trace particles reduces to do (left); the connection of 
the particles on a magnetic flux tube changes after the reconnection 
(right). 



A. Shear flow test 




FIG. 6: Snapshots of two simultaneous reconnection events, before 
the reconnection (top) and after (bottom). Note the change of con- 
nections of the flux ropes after the reconnection. 



We find excellent agreement between Eq. © and our numeri- 
cal solution. In Fig. [3] we plot |B| at a fixed value of y versus 
time to compare it with Eq. ©. The comparison is quite satis- 
factory; the step-wise change in the numerical solution for B 
arises because, in this simple flow with a the shear rate slowly 
varying in space, many new trace particles are introduced si- 
multaneously as the flux tube is stretched, and then no parti- 
cles are added for some time until the next series of particle 
insertions. After t w 20, we reverse the flow to observe that 
the particles and magnetic field return to their initial states, 
to demonstrate that our algorithm correctly describes the con- 
traction of the flux tubes as well. Fig. [4] shows that the flux 
tube adopts the shape of the flow before the flow field is re- 
versed, colour coding indicating the magnetic field strength. 



In order to test our model we consider a two dimensional 
shear flow with a Gaussian profile, 

u=(u x ,0), u x = u e~ v2/2 , (3) 

and a flux tube extended across the flow from y = — oo to 
+oo. For £B = 0, Eq. (fl} can easily be solved exactly to 
yield, 

|B| = B ^l + u y 2 e-y 2 t 2 , (4) 
where Bq is the initial field strength. Since, 



\B\dV = 



ip(l) dl cx L, 



(5) 



where ip = const is the magnetic flux and L is the length of 
the flux tube, and since B is independent of x, we have 



LocB Vl + u y 2 e-y 2 t' 2 dy. 



(6) 



B. Reconnections 

Reconnections are introduced into the model in a straight- 
forward manner. If the separation between two particles, 
which are not neighbours, becomes less than a certain scale 
do, we reconnect their associated flux tubes by reassigning the 
flags P which identify the particles ahead and behind those in- 
volved in the reconnection, as shown in Fig. [5] We found that 
do has to be comparable to the separation of the trace parti- 
cles, d, in order to obtain meaningful numerical results, e.g., 
do = 1.5c?. Two particles are removed from the system after 
each reconnection event, those labelled P = 2 and P = 12 
in Fig. [5] and their magnetic energy is lost, presumably to 
heat. We also monitor the cross product of the tube tangent 
vectors close to the reconnection point. By ensuring that the 
magnitude of the cross product is smaller than some tolerance 
e w 10~ 2 and that the magnetic fields in the reconnecting 
loops are (almost) oppositely directed, we prevent parallel flux 



4 



» 





t=5.0 























-2 


-1 1 

X 


2 




t= 


40.0 
























FIG. 7: (Colour online) Flux expulsion by differential rotation Q: 
the form of magnetic line initially aligned with the y-axis at different 
times indicated in the corner of each frame. The field strength grows 
as the magnetic line is wound around by the differential rotation. 
Eventually the separation of neighbouring turns becomes less than 
do and reconnections destroy the field. Magnetic field strength is 
colour coded as in Fig. [4] 




150 



I 4.4 



3.8 
3.6 



4 - 

log R„ 



FIG. 8: The upper panel shows the root-mean-square magnetic field 
strength B !ms as a function of time for the simulation shown in Fig. [7] 
The lower panel represents the scaling of the maximum values of 
firms among eight simulations with decreasing do. The line of best 
fit, shown dotted, has the slope 0.54 ± 0.21 



tubes from reconnecting. In Fig. [6] we show snapshots from 
a simulation before and after two simultaneous reconnection 
events. Since we monitor the amount of magnetic energy re- 
leased in each reconnection event, we know the total magnetic 
energy released by the reconnections over any given time pe- 
riod. Finally we introduce a minimum loop size of 3d, i.e., 
no magnetic loop can contain less than three particles. Any 
smaller loop is removed from the system, releasing its energy. 
We shall see later that this cutoff is important when we take 
derivatives along the loops to calculate magnetic tension. 



C. Flux expulsion 

We test the reconnection algorithm by considering mag- 
netic flux expulsion from a region with closed streamlines 
lfl7ll . Consider an initially uniform magnetic field Bo, in our 
case a single flux tube extended over — oo < y < oo along 
x = 0. Differential rotation is applied to the field, with veloc- 
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ity given in cylindrical polar coordinates by 
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U = (u r , Ug), U r = 0, U0 
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Solutions of the induction equation grow linearly in time until 
a maximum magnetic field is achieved, 



|B| max = O« 2 )B , 



(8) 



where R m is the magnetic Reynolds number, after which mag- 
netic diffusion destroys the field in the rotating region. We 
find a similar scaling of the maximum magnetic field strength, 



with the dimensionless quantity 



Rm — 



u r do ' 



(9) 



which we identify as the effective magnetic Reynolds number. 
Here uq and Iq are typical velocity and length scales respec- 
tively, d is the reconnection length, and u r is the character- 
istic reconnection speed. In the case of the Gaussian vortex 
this is taken as the relative velocity of the approaching flux 
tube. Figure [7] shows snapshots of a typical simulation as it 
proceeds: the magnetic field after one winding (t = 5.0), 
in a state close to the maximum field strength {t = 40.0), 
as the reconnections start to drive the destruction of the field 
(t = 46.5), and finally the quasi-steady state (t = 200.0). The 
first plot in Fig.|8]shows the corresponding values of B Ims ver- 
sus time; the linear growth before the onset of reconnections 
is apparent. The second plot in Fig. [8] shows the power-law 
relationship between i? rms ,max and R m : the slope of the fit 
shown is 0.54 ± 0.21, in a reasonable agreement with Eq. ([8]). 
One final test, results not presented here, was to ensure that 
no dynamo could be supported by driving the flux ropes with 
a two dimensional flow, i. e. u = (u x , u y ,0). 



III. MODEL OF A TURBULENT FLOW 

Our next step is to choose the velocity field which drives the 
dynamo. Following previous work [28], to bypass the com- 
putational limitations of direct numerical simulations (DNS), 
we use the so-called Kinematic Simulation (KS) model. The 
KS model has primarily been used as a Lagrangian model 
of turbulence and results are in good agreement with DNS 
ami 11811 . This flow is known to be a hydromagnetic dy- 
namo [28]. The KS model prescribes the flow velocity at a 
position x and time t through the summation of Fourier modes 
with randomly chosen parameters (which are then kept fixed), 



according to 
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u(x, t) = (A n x k„ cos 4>n + B„ x k n sin < 



(10) 



where <j) n = k„ • x + u> n t, N is the number of modes, k ra and 
w n = k n u n are their wave vectors and frequencies (see ll28"ll 
for details). The unit vectors k n are chosen randomly, and 
k n = fc n k n where k n is the wave number of the n'th mode. 



We choose the directions of A„, and B„ randomly, imposing 
only orthogonality with k„ , which gives 



A.ri x k« 
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and likewise for B„. We then select a kinetic energy spectrum 
E(k) and set 



An = B n = [§£(fc„)Afc„] 
This ensures that 



1/2 



(12) 



i p pOC _____ 

- ±\u\ 2 dV = E{k)dk^Y^E{k n )Ak n . (13) 

V Jv J ° n=l 

One of the main advantages of the KS model is that we have 
complete control of the energy spectrum, E{k n ) via appro- 
priate choice of A„ and B„. We also note that V • u = 
by construction. We adopt a modification of the von Karman 
spectrum, 



E(k) = k\l + fc 2 ) -(2 +P /2) e -i( fc/fcN ) 2 



(14) 



which reduces to E(k) oc k~ p in the inertial range, 1 <C 
k -C fcjv, with k = 1 at the integral scale; p = 5/3 pro- 
duces the Kolmogorov spectrum, and k^ is the cut-off scale. 
Figure [9] shows the energy spectrum of the KS flow, obtained 
numerically after fast Fourier transforming u calculated from 
Eq. ( fTOb on a 128 3 mesh. We also show a slice, in the (x, y)- 
plane, of the corresponding vorticity field. 

The results presented below have been obtained with k\ = 
2ir and k^ = 167T, so that the smallest velocity scale is 
In = 2ir/kN — 0.125. For comparison, the reconnection 
scale is adopted as do = unless stated otherwise. With 
this prescription, the effective magnetic Prandtl number in our 
model is larger than unity. 

To check if our results depend on the form of the flow, we 
also use an ABC flow of the form [7], 

u= (cosky+sinkz,smkx+coskz,coskx+smky) , (15) 

known as the 111 ABC flow. Dynamo action driven by ABC 
flows has been studied extensively [10]. This particular flow 
has eight stagnation points for < k x < 2ir. In the plane 
in which the flow converges to a particular stagnation point, 
the magnetic field is advected, and becomes elongated in the 
direction of the streamlines which diverge from the stagnation 
point. The resulting magnetic structures are commonly de- 
scribed as 'magnetic flux cigars' |H. Figure [10] shows such 
magnetic structures produced by our numerical solution of 
the induction equation, along with the corresponding veloc- 
ity field. 



IV. DIFFUSIVE AND RECONNECTION-BASED 
DYNAMOS 

Comparisons of the solution of the induction equation with 
those produced by the flux rope model are not straightforward 
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FIG. 9: (a) The energy spectrum E(k) as obtained by Fourier transform of Eq. l llOt with N = 20, ki — 10 and fcjv = 400. The dashed line 
has E(k) oc fc _5//3 . (b) Slice in the (x, j/)-plane of the vorticity field from (a), lighter shades indicating higher vorticity. Velocity vectors are 
shown in white. 




FIG. 10: (Colour online) The distribution of |u| in the 111 ABC 
flow Eq. H5\ is shown in colour, with the magnetic isosurfaces with 
|B| = 3.5-Brms, obtained by solving the induction equation, shown 
in grey scale. Each point on a 3D mesh is assigned an opacity and 
colour, depending on |u|. Regions with |u| > 2.5urms are coloured 
purple and those where |u| ft; are coloured blue. 



because of the difference in the control parameters of the two 
models: the magnetic Reynolds number R m = uqIq/t) and 
the reconnection length do, respectively. Following the pre- 
vious section, we introduce the effective magnetic Reynolds 
number as R m = uqIq/ (iirdo), where u r is the characteristic 
reconnection speed. 

We find that the dynamo based on reconnections is more 
efficient than the diffusion-based dynamo, in the sense that 
the growth rate of magnetic field of the former is significantly 



larger when R m w R m . Therefore, in order to achieve conser- 
vative conclusions, we compare dynamos with similar growth 
rates of magnetic field. Thus, we have R m > R m in the mod- 
els which we compare. 

An advantageous property of both KS and ABC flows 
is their analytic nature, which means that we can follow 
(Lagrange-like) fluid particles in the flow without using an 
Eularian mesh. The initial condition of our simulations is an 
ensemble of random closed magnetic loops; both the induc- 
tion equation and the flux rope model are evolved with the 
same velocity field (apart from the overall normalisation to 
provide comparable growth rates of magnetic field). The ini- 
tial condition for the induction equation is obtained by Gaus- 
sian smoothing of the magnetic field in the ropes, where we 
define the smoothed field B(x), as 

6(x) = 2^J V eHx ~ y|2/2CT2]3 (y) d y 3 - < 16 ) 

Importantly, this procedure preserves the solenoidality of the 
field, i. e., V • B = 0. Figure [TT| shows the smoothed ini- 
tial magnetic field used in a typical simulation, along with the 
corresponding flux rope setup. 

The induction equation is solved using the Pencil Code JUl, 
which implements a high-order finite-difference scheme, on a 
256 3 mesh with 1000 < R m < 1500 in a periodic box. Sim- 
ulations with the KS velocity field had fci = 2ir, Un — 1Q"K, 
and p = —5/3; here i? m is based on the largest velocity scale 
2n/ki. In a separate simulation flux ropes are advected and 
stretched by the same velocity field, where the positions of 
the trace particles are evolved using a 4 th order Runge-Kutta 
scheme, with a time step of In / (20u at). The algorithm for in- 
serting and removing particles is applied every time step, and 
the reconnection algorithm, every ten time steps. We choose 
d to be 1/4 of the smallest length scale in the flow and set 
do/d — 1.5, where g?q is the reconnection length scale. 
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FIG. 1 1 : (Colour online) Initial conditions for the magnetic field in 
the ABC and KS simulations. The upper panel shows the initial set of 
flux ropes populated by the trace particles. The lower panel presents 
a direct volume rendering of a similar set of closed flux ropes, but 
now smoothed using a Gaussian kernel d!6l l. The smoothed magnetic 
field is used to initialise the induction equation simulations. 



In Fig. Q~2] we present snapshots of the magnetic field from 
the flux rope simulations, driven by the KS flow, as it evolves, 
Fig. [13] shows the corresponding evolution in the ABC flow. 
In the case of the ABC flow, one notices an anticorrelation 
between the curvature of the flux tube and the magnetic field 
strength (colour coded) as suggested by Schekochihin et al. 
JH. 

Our model is deliberately oversimplified with respect to the 
(incompletely understood) physics of magnetic reconnections. 
Nevertheless, we can argue that our model is conservative 
with respect to the reconnection efficiency. The reconnecting 
segments of magnetic lines in our model approach each other 
at a speed u r ~ uoRe -1 ' 4 for the Kolmogorov spectrum, 
equal to velocity at the small scale do <C lo with lo the energy- 
range scale of the flow and do assumed to be close to the tur- 
bulent cut-off scale. If the magnetic field is strong enough, 
the Alfven speed Va, which controls magnetic reconnection 
in more realistic models, is of order u(Iq). Then u x <C Va 



t=0.09 




FIG. 12: (Colour online) The evolution of magnetic flux tubes in the 
KS flow: snapshots taken at times given in the upper right corner of 
each panel. Magnetic field strength is colour coded, with the colour 
bar shown next to the last snapshot. Note the overall increase of 
magnetic field strength as time proceeds. 
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FIG. 13: (Colour online) As in Fig. [12] but for the 111 ABC flow. 
Magnetic field has been stretched into 'flux cigars' which are even 
more apparent when the field is later smoothed. 



and our model is likely to underestimate the efficiency of re- 
connections. The Sweet-Parker reconnection proceeds at a 

— 1/2 

speed of order Va-Rui . whereas the Petschek reconnection 
speed is comparable to VA/lni? m [20]. For uq ~ V\ and 
i? m ~ Re ^> 1, the reconnection rate in our model is larger 
than the former but much smaller than the latter. 



COUPLING WITH THE NAVIER-STOKES EQUATION 
AND DYNAMO SATURATION 



A remarkable feature of the model of reconnecting mag- 
netic flux tubes suggested here is that is admits straightfor- 
ward extension to include the back-reaction of magnetic field 
on the flow via the Lorentz force. Such a generalisation is pre- 
sented in this section. To obtain a flow similar to the KS flow 
( fTOb as a solution of the Navier-Stokes equation, we include a 
driving force shown here after the viscous term: 



Du 
~~Dt 



VP 

P 



J x B 

47T/9 



1 

Re 



V 2 u 



u K s 



(17) 



where D/Dt = d/dt + u ■ V is the convective (Lagrangian) 
derivative, Re is the Reynolds number, uks is the KS velocity 
field ( [Tol l, and r is a certain relaxation time. The smaller is r, 
the closer u is to the KS flow. Since we assume that magnetic 
field is localised within flux ropes, magnetic pressure must be 
balanced by some other force, presumably by gas pressure, so 



we assume that V(P + B 



0, and only the magnetic 



tension force (B • V)B remains to be balanced in the Navier- 
Stokes equation. Neglecting viscosity, Re — > oo, we then 
obtain 



^ = -(B-V)B 

Dt 8tt v ' 



uks 



(18) 



If magnetic field is confined into thin ropes and aligned with 
their axes, magnetic tension involves the directional deriva- 
tive of magnetic field along the rope axis alone, (B • V)B = 
B 9B / ds, where s is the distance measured along the rope. 
Thus, it is sufficient to have magnetic field defined on mag- 
netic loops (rather than at any position in the volume) in order 
to calculate magnetic tension force. 

We require a solution of Eq. ( [T8l at the changing positions 
of the trace particles, so we need, essentially, a Lagrangian so- 
lution of this equation. Assuming that the flow is close to the 
relaxed state and does not change rapidly, we put Du/Dt ss 
to obtain the trace particle velocities as 



u 



<9B 

u K s + tB—. 

OS 



(19) 



This approximation filters out rapid wave motions, e.g., Al- 
fven waves, which simplifies numerical simulations. We use 
this approximation to study the saturation of the dynamo ac- 
tion in Section lVAl where we consider rather long time inter- 
vals. 

On the other hand, our model also allows us to include Al- 
fven waves and their nonlinear interactions. For this purpose 
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we assume that \Du/Dt\ |uks — u|/t and the Navier- 
Stokes (or rather Euler) equation reduces to 



Dm 

~Dt 



B 



.OB 

ds ' 



(20) 



which, can be coupled with the equation for a frozen-in mag- 
netic field DB/Dt = (B • V)u, written in a similar form: 



dB 

~dt 



OS 



(21) 



Imposing a homogeneous magnetic field Bo and a weak per- 
turbation, and linearising these equations leads to the wave 
equation describing the Alfven waves. Since we assume that 
magnetic pressure is precisely balanced by gas pressure, our 
model does not admit compressible waves. 

The nonlinearity requires that we make two changes to our 
numerical calculations. Firstly in our model B dB / ds is only 
defined at positions on the magnetic line (flux tube), and so 
the velocity field can only be evolved at those positions. The 
fourth-order Runge-Kutta time stepping scheme used in the 
kinematic regime is not suitable as it requires velocity field at 
positions where magnetic field is not defined. Therefore, we 
use the three-step Adams-Bashforth scheme instead to evolve 
the positions of the trace particles: 

h 

x„ +1 = x„ + — (23u„ - 16u„_i + 5u n _ 2 ), 

t n +l = t + h, 

where h is the size of the timestep. 

The differentiation of the magnetic field along the flux tubes 
requires an improved accuracy for the positions of newly in- 
troduced trace particles in a stretched flux tube. A first-order 
prescription (O is no longer accurate enough and we replace 
it by a second-order interpolation scheme. Consider a section 
of magnetic line traced by three particles at positions xi, x 2 
and X3. If the distance between x 2 and X3 becomes greater 
than d, we introduce a new particle at the position X4 given by 

x 4 = xi - [(x 3 - xi) - 4(x 2 - xi)]/x 
+ [2(x3 - xi) -4(x 2 -Xi)]^ 2 , 

where /j, is a parameter. For /j, = 0.75, the new particle is 
placed between x 2 and X3 as required. In tests, in particu- 
lar with Alfven waves, we found a substantial improvement 
in the accuracy of the solution with this higher-order scheme. 
However kinematic results show no quantifiable difference be- 
tween the two schemes. 

The directional derivative of magnetic field has to be calcu- 
lated carefully in our case since the separation of trace parti- 
cles is not constant. We use the following numerical schemes 
to evaluate the first derivative: 
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0(f), (22) 



where the notation is defined in Fig. [T4J with r replacing B. 
For li = ti_i = h, we recover a commonly used finite differ- 
ence scheme. 
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FIG. 14: A sketch of tracer particles that populate the flux tubes 
illustrating Eqs. {22} and d24D 
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FIG. 15: The evolution of the root-mean-square (r.m.s.) magnetic 
field B Ilns in nonlinear dynamo driven by the KS flow d 1 Ob with 
the nonlinearity (19}. Exponential growth at t < 4 is followed by 
a saturated state where magnetic energy density fluctuates around a 
roughly constant level. The unit time is the kinematic time scale at 
the largest scale in the flow, to = Io/uq. 



A. Saturated dynamos 

Our starting point here is Eq. ( fl9] l for the velocity field. At 
each position on a flux tube, x", we calculate the KS veloc- 
ity field using Eq. ([Tol l, and modify it with magnetic tension 
force. The directional derivative of magnetic field along the 
tube, dB/ds, is computed using Eq. (l22l . 

The details of the simulations are similar to those in the 
kinematic regime. We find our timestep of Zjv/(20mat) to be 
sufficient to capture the dynamics of the motion (we found no 
noticeable difference between simulations with the timestep 
set an order of magnitude smaller than this). The relaxation 
time t is set to be the same as the timestep, r = O(10~ 5 ). 

Figure [15] shows the r.m.s. magnetic field strength B lms as 
a function of time, where the initial exponential growth is fol- 
lowed, at t > 4, by the saturation of the dynamo action, with 
B rms fluctuating around a roughly constant level. 

Snapshots of the magnetic ropes shown in Fig.[T6lillustrate 
the spatial evolution of the magnetised region. Unlike vir- 
tually all earlier simulations of the fluctuation dynamo, most 
often performed in periodic boxes with volume-filling initial 
conditions, the initial magnetic field in our simulations is lo- 
calised in space, as shown in the upper left panel of Fig. [16] 
(which refers to an early stage of evolution). In the kine- 
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FIG. 16: (Colour online) Snapshots of the magnetic loops evolving in the KS flow l llOb in both the kinematic stage t < 4 and in the saturated 
state t > 4, taken from the same run as Fig.[T5] Magnetic field strength in the flux tubes is colour coded, with the colour scheme shown in the 
bottom right corner. Note the scale of the box remains the same (dimensions shown at t = 0.6) for all snapshots. Notice the high density of 
the snapshot when t = 3.9 and the field is at a maximum. This corresponds to an overshoot in B rms visible in Fig.[l5]as the dynamo saturates. 



matic stage, t < 4, magnetic field growth is accompanied by 
the spread of the magnetised region clearly visible in the first 
three snapshots. Consistently with the action of magnetic dif- 
fusion ?7t oc IqUo, the size of the region occupied by magnetic 
ropes grows at t 1 / 2 However, the spread is halted in the non- 
linear, saturated stage represented in Fig. [16] by snapshots at 



t = 25, 50 and 85, which suggests that the turbulent mag- 
netic diffusivity is suppressed in the saturated state. This ap- 
pears to be a result of the suppression, by magnetic tension, 
of random stretching of magnetic field at the location of a 
flux rope. In other words, the saturation of the dynamo ac- 
tion is achieved via the suppression of the effective magnetic 
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FIG. 17: As, the average separation of trace particles arranged along 
magnetic loops in the KS flow d 10b . grows rapidly in the kinematic 
regime of the dynamo, and then remains roughly constant as soon 
as the dynamo action saturates at t > 4 when the magnetic tension 
force becomes significant. 



Reynolds number, R m = IqUq/t), where r) is the effective mi- 
croscopic magnetic diffusivity. This idea is fully consistent 
with the arguments of Subramanian B26I1 who considered a 
similar nonlinearity in the Kazantsev model of the fluctuation 
dynamo [see also [l3l [Til f27ll and also of Schekochihin et al. 
B23l 12411 who suggested a simple model of the effect of Lorentz 
force. A general feature of these models is that the small-scale 
dynamo saturates because of a 'renormalization' of the co- 
efficients governing its evolution, and the corresponding de- 
crease in the effective magnetic Reynolds number. This can 
be the result of enhanced nonlinear diffusion B26I1 . increased 
diffusion together with additional hyperdiffusion If27ll . or re- 
duced stretching |[T3i[T4ll23il24ll . Our model is consistent with 
the saturation of the dynamo action via the suppression of the 
magnetic Reynolds number, now arising from a reduction of 
localised random stretching, or turbulent magnetic diffusivity. 

To clarify further the mechanism of dynamo saturation, we 
monitored the average separation of trace particles in the flow. 
At the start of a simulation each trace particle located on a 
magnetic loop is assigned a neighbour, for convenience the 
particle next to it. As the simulation proceeds, the particles 
are advected by the flow, and new particles may be introduced 
between them, but we continue to monitor the separation be- 
tween the original pair of particles. The particle separation 
averaged over all the original particle pairs, A#, is shown in 
Fig. [17] Indeed, the separation of the particles stops growing 
as soon as the dynamo enters the nonlinear stage. We stress 
that, at late times, not all pairs of trace particles belong to the 
same magnetic loop because of multiple reconnections that 
often split a magnetic loop into smaller ones. Thus, the fact 
that A b ceases to grow implies that not only stretching is sup- 
pressed within a single loop, but also that the magnetic loops 
stop spreading in space. 

We performed another experiment, where the location of 
the test particles used to compute the dispersive properties of 



3 




FIG. 18: A„, the average separation of 1000 test particles in snap- 
shots of the KS flow JlOb . unaffected (solid) and affected (dashed) by 
magnetic field. 



the flow was not restricted to the magnetic loops. To reduce 
technical problems, we considered two time-independent 
flows obtained as the snapshots of the original KS velocity 
field and of its form affected by magnetic tension force at 
a certain moment in the saturated dynamo state. The evo- 
lution of the particle separation (averaged over 500 pairs) is 
shown in Fig. [18] where one can distinguish the initial expo- 
nential growth of the separation, followed by the Richardson 
regime A n oc i 3 / 2 when A u < l Q , and then by the incoherent, 
diffusive dispersion A u oc t x l 2 at larger separations. How- 
ever, the separation of the particles is insensitive to the effects 
of the Lorentz force: the two curves in Fig. [18] hardly differ 
from each other. Thus, the flow has not been strongly affected 
by the Lorentz force, except for the close vicinity of the flux 
tubes. 

Importantly, the mean particle separation A u is a measure 
of kinetic turbulent diffusivity, as opposed to the magnetic 
one. The former involves the mean square of the total ve- 
locity v t oc (tu 2 ), whereas the turbulent magnetic diffusiv- 
ity is only sensitive to the velocity field components orthog- 
onal to the magnetic field, v% oc (ritj_). Incidentally, these 
results imply that the turbulent magnetic Prandtl number is 
different from unity: in isotropic flow and magnetic field, 
Pr m = v t /rjt ~ 3/2; the difference is small but perhaps sig- 
nificant in some applications. We cannot exclude that this fea- 
ture is an artifact of our model where a localised modification 
of the velocity field by magnetic tension does not spread into 
a broader region as it would do due to kinematic viscosity. 

B. Alfven waves 

To demonstrate the flexibility of the flux model suggested 
here, we briefly present simulations of Alfven waves prop- 
agating along an imposed magnetic field. For this purpose 
we solve Eqs. (l20l i and d2ll with a single magnetic flux tube 
which, initially, has constant magnetic field. The tube is per- 
turbed as shown in the upper left panel of Fig. [19] The sim- 
ulation proceeds in a box periodic in the x-direction, and the 
interaction with the ghost wave from the next periodicity cell 
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FIG. 19: Snapshots of the Alfven wave simulations. A single flux tube with field strength |B| = 1 is perturbed at the midpoint. As the 
simulation progresses two wavefronts form which move apart with speed equal to the field strength. As the fronts reach the edge of the 
periodic box they interact at the boundary. 



is seen in the final snapshot at t = 14.8. We confirmed that the 
phase speed of the wave is indeed proportional to the strength 
of the magnetic field. Since our model admits nonlinear in- 
teractions of Alfven waves, it can be used to study spectral 
energy cascades and other features of the Alfven wave turbu- 
lence. 



VI. CURVATURE OF MAGNETIC LINES 

Schekochihin et al. [22] discuss the geometry of magnetic 
lines in the kinematic fluctuation dynamo driven by a single- 
scale, (5-correlated in time random flow with high magnetic 
Prandtl number. They argue that magnetic field strength and 
magnetic line curvature should be anticorrelated and derive 
a power-law probability distribution function of the magnetic 
line curvature. These results are used to support the picture 
of folded magnetic lines as a representation of magnetic field 
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FIG. 20: The scatter plot of curvature, k versus |B| calculated at the 
end of a simulation, with both variables normalised to the respective 
maximum values. 
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FIG. 21: Probability density functions (PDF) of curvature, k at a late 
stage of magnetic field evolution in the KS flow illustrated in Fig.l 121 
Dashed line is for the power-law distribution P(k) oc k _13//t ob- 
tained by Schekochihin et al. 12211 . 



produced by the fluctuation dynamo. 

The anticorrelation between the curvature of magnetic lines 
and the strength of the magnetic field is intuitively appealing 
since magnetic field strength grows due to a random stretch- 
ing of magnetic lines which is necessarily accompanied by 
a reduction in their local curvature. However, the stretching 
is not the only component of the fluctuation dynamo mecha- 
nism. In the framework of the stretch-twist-fold dynamo con- 
cept, stretching must be followed by the folding of magnetic 
lines to ensure an exponential growth of magnetic field - and 
the folding will tend to increase the local magnetic line cur- 
vature. Therefore, an anticorrrelation between magnetic cur- 
vature and strength may be expected for a decaying magnetic 
field rather than for magnetic fields growing due to the dy- 
namo action. In this section we explore directly the relation 
between the magnetic line curvature and strength using the 
reconnecting flux rope dynamo model. 

The curvature of the flux ropes can be calculated as Hill 



(23) 
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where r(s) is a parametrised space curve representing a mag- 
netic flux rope, with s the distance measured along the rope, 
and dash denotes derivative with respect to s, with the first 
derivative calculated using Eqs. d22b and the second deriva- 
tive, from 
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2r, ; 



0{f), 
(24) 

where notation is defined in Fig. [14] For = = h, this 
form reduces to a standard finite difference scheme. 

Figure[20]is the scatter plot of magnetic field strength at a 
particular position versus the field line curvature at that point, 
computed at the end of the simulation, illustrated in Fig. [12J 
where the dynamo is driven by the KS flow. Only the envelope 
of the curvature distribution appears to be consistent with the 



anticorrelation, and even that only for relatively strong fields, 
B > 10~ 2 B max . Thus, the range of the curvature values is 
narrower at positions where the field is stronger, but for any 
field strength this range includes very small curvature values. 

Those parts of magnetic flux tubes where magnetic field is 
weak have low curvature, especially those with B <C Bo with 
B the initial field strength (B w 4 x 10~ 3 -B max at the par- 
ticular time of the simulation). In our model, the only way 
the field strength |B| can become smaller than Bo is through 
the shrinking of a flux tube caused by contracting flow. In 
a perfectly conducting fluid, such a contraction can make the 
curvature larger, e.g., when a wavy magnetic line is contracted 
along its wave vector. However, the situation changes entirely 
in the presence of reconnections (or any other magnetic dissi- 
pation mechanism): now, reconnections eventually eliminate 
the bends of the magnetic line thus reducing the curvature of 
a contracting magnetic line. Apparently, we see the evidence 
of this in Fig. [20] The group of points with nearly maximum 
curvature at the top of the frame are probably those which will 
undergo reconnections of this type very soon. Finally we note 
that the reconnection length do limits the maximum value that 
k can take. 

Schekochihin et al. [22] showed, both analytically and nu- 
merically, that the probability density function P of the cur- 
vature of field lines has a power-law form in the limit of large 
k. In particular, they obtain P(n) oc k~ 13 / 7 for a three- 
dimensional, incompressible flow. Figure |2"TI shows the PDF 
of curvature from our simulations which shows a very good 
agreement with the analytical results of Schekochihin et al. 
& 

Schekochihin et al. [25] argue that, because the magnetic 
field strength is higher where the magnetic line curvature is 
lower (due to the stretching by velocity shear), magnetic ten- 
sion force (which is quadratic in B) is lower at positions with 
weak field. This conclusion relies on estimating the tension 
force as \BdTS/ds\ ~ kB 2 , and noting that this quantity 
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FIG. 22: The slope of the PDF of magnetic line curvature and 
r. m. s. magnetic field in the kinematic and saturated states of the 
dynamo. The solid line shows a, where PDF(k) ~ k°. In the kine- 
matic regime a w —13/7 as in Fig. [2T| As the dynamo saturates 
(log B rms shown dashed) a decreases. 



decreases with B 2 , although n increases. This would imply 
that magnetic tension is relatively unimportant in regions with 
strong curvature. This leads these authors to a conclusion that 
magnetic field is organised into folded structures which persist 
in the saturated state. However, what matters is not the value 
of B 2 itself but rather its gradient along the magnetic line, 
dB 2 /ds. In a stretched magnetic line, the gradient is reduced 
in regions with strong magnetic field (i.e., the straight seg- 
ments of the folded magnetic lines) and enhanced in regions of 
weaker field (in the turns in the folded structures). Therefore, 
magnetic tension will drive the turns closer to each other along 
each magnetic line destroying the folded structures. Thus, the 
lack of any apparent domination of folded structures in Fig. [12] 
is consistent with the curvature PDF shown in Fig. [21] We 
show in Fig. [22] the time variation of the slope of the curva- 
ture PDF into the nonlinear regime (discussed in Section [V]): 
the PDF becomes steeper , so that high curvature occurs less 
often in the nonlinear state. This can be attributed to magnetic 
tension which tends to reduce magnetic line curvature. 



VII. STATISTICS OF MAGNETIC ENERGY RELEASE 

Solar corona is one of the astrophysical environments 
where magnetic reconnections are believed to play important 
role, particularly in heating the plasma to the high temper- 
atures observed [20]. The reconnections are assumed to be 
driven by the motion of the footpoints of magnetic flux tubes 
anchored in the photosphere and extending into the corona 
Ir2ll and references therein]. Reconnection events that release 
large amounts of magnetic energy are observed as solar flares. 
A remarkable feature of the coronal heating mechanism is that 
the frequency distribution of the flare energy has a power law 
form in a very broad energy range (eight orders of magnitude) 
[see an excellent review of|6J] 



P(AM) oc (AM)' 



If s < —2, most of the magnetic energy released into the 
corona is due to weak flares. This attractive option suggested 
by Parker lfl9ll is known as the nanoflare model of the coronal 
heating. This idea is most often explored in the context of self- 
organised criticality models based on cellular automata, which 
are known to demonstrate the required power-law statistical 
distributions. Notably, the continuous analogies of these mod- 
els involve the hyperdiffusion operator [6]. A widely recog- 
nised difficulty of this approach is the elusive connection with 
the physical picture and even unclear physical interpretation 
of the variables. Alternative models [e.g., [HI], where recon- 
nection evens are modelled directly, also reproduce the power- 
law statistics, but still remain rather idealised regarding the 
behaviour of magnetic flux tubes. 

Our model is quite different from the Solar corona settings, 
where the reconnections are driven by the motion of the flux 
rope footpoints, the plasma is believed to be magnetically 
dominated, and in situ dynamo action is improbable. Never- 
theless, in this section we consider the statistics of the energy 
release in our model of the flux rope dynamo. As we show 
here, our reconnection dynamo model naturally develops a 
power-law distribution ( FSB! with s ~ — 3, which appears to be 
independent of the form of the velocity field. Our model can 
readily be adapted to the Solar corona conditions, and despite 
the differences of our model from the Solar corona models, 
we feel that this feature of the model can be relevant in this 
context. 

In the case of the induction equation, the magnetic energy 
dissipation rate can be defined as 
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where M is the total magnetic energy. A similar quantity can 
be obtained for the flux rope dynamo by summing the con- 
tributions of all reconnection events to the magnetic energy 
release: 
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where r is a suitable time interval during which N T reconnec- 
tions occur (we take r to be equal to ten time steps; individ- 
ual reconnection events occur in a single time step), and B i7 
Si and Li are the magnetic field strength, the cross-sectional 
area and length of the reconnected (and thus removed) flux 
tube segment associated with a trace particle number i. From 
our assumption of frozen flux BiSi = ip — const, the total 
magnetic energy M is, 



B 



(28) 



where A^ot is the total number of trace particles in all flux 
tubes. Thus, 
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FIG. 23: Magnetic energy release rates from two kinematic dynamo 
models with the KS flow and similar growth rates of magnetic field: 
as obtained from the induction equation (black) and the reconnecting 
flux rope model (grey). The former has a mean value of 2.4 (here 
R m = 1200) once the eigensolution has developed. The latter (with 
R m — 174) has a mean value of 23 (shown with thick white horizon- 
tal line). 
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FIG. 25: Probability density for the scaled magnetic energy release, 
£ = (7 — 7")/cr T , from the time series of Fig. [23] for the flux rope 
dynamo (circles) and the diffusive dynamo with the same magnetic 
field growth rate and velocity field of the same form (squares). A 
power-law fit to the former and a Gaussian fit to the latter are shown 
solid and dashed, respectively. 
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FIG. 24: As in Fig. [23] but for the ABC flow with ydo/u « 0.6 
for the eigensolution of the induction equation at R m = 55, and 
7r?o/wo ~ 6.7 with R m — 24 in the flux rope dynamo. 
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FIG. 26: As in Fig. [25] but from the time series of Fig. [24] a log- 
log plot for the flux rope dynamo (circles) with solid line having the 
slope —2.98. As above, we show a Gaussian fit (dashed) to the data 
from the diffusive dynamo (stars) driven by the ABC flow. 



In Fig.[23]we present the energy release rates in simulations 
where the growth rate of the magnetic field is a = 0.16inboth 
simulations (with the unit time Iq/uq). The dashed line shows 
the energy release rate from a simulation of induction equa- 
tion with R m = 1200, which has the mean energy release rate 
7i 2.4. The solid line shows the corresponding results from 
the reconnection dynamo, with the mean value plotted as a 
thick horizontal line. The mean value of the energy release 
rate from the reconnecting flux rope dynamo is j r w 23, an 
order of magnitude larger than that obtained from the induc- 
tion equation. We also note the strong fluctuations in the en- 
ergy release rate from the reconnection model, as opposed to 
the quiescent behaviour in the induction equation. It is impor- 
tant that an order of magnitude difference in the energy release 



rates occurs in solutions with similar growth rates of magnetic 
field. Since the reconnection based dynamo is more efficient 
than that based on magnetic diffusion (see section HVb kinetic 
energy density in the former being 10 times smaller than in 
the latter. With comparable kinetic energy densities, the dif- 
ference between the energy release rates can be even larger. 

As shown in Fig. [24] dynamos driven by the ABC flow be- 
have similarly. With R m = 55, the induction equation gives 
an energy release rate of about 71 = 0.6. The corresponding 
flux rope dynamo with the same growth rate (er = 0.02) has 
energy release rate 7; w 6.7, again ten times larger. 

We show in Figs. |25] and |26l the probability distributions of 
the magnetic energy release rate, normalised to the total mag- 
netic energy in the domain, £ = (7 — 7)/<x 7 , with 7 = 7; 
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FIG. 27: As in Fig. [25] but for the energy release rates from a non- 
linear simulation (see Fig.ll5t. A power law fit is shown with dashed 
line (with slope —3.1) 



or 7 r . Here overbar denotes time averaging (at times where 
an eigensolution has been established) and 07 is the standard 
deviation of 7. Since 7 = AM / (Mr), it can easily be seen 
that C = (AM - AM) /a am- We obtained the probability 
distributions of £ from both the induction equation and the re- 
connection dynamo model, both driven by the KS flow, shown 
in Fig. [25] and also with both based on the 111 ABC flow, 
shown in Fig. [26] The power-law index obtained for the KS 
flow is s « -3.3, and that for the ABC flow is s ps -3.0 . We 
stress that this power-law behaviour is not related to the na- 
ture of the velocity field: solutions of the induction equation 
with the same velocity fields, as we show in Fig. [25] exhibit 
an approximately Gaussian probability distribution. 

Results shown in Figs. [25] and [26] have been obtained from 
kinematic simulations, where the velocity field was not af- 
fected by the Lorentz force. However, the corresponding non- 
linear model introduced in Section[V]retains this feature, with 
s w —3.1 for the KS flow in the statistically steady state, as 
we show in Fig. [27] 

It is not quite clear if the flux rope dynamo represents a 
physical example of self-organised criticality, but the system 
seems to possess at least some of the required properties. In 
particular, as we argue above, our model can be viewed as an 
extreme case of magnetic hyperdiffusivity, which also arises 
in the self-organised criticality models of Solar flares. 



VIII. DISCUSSION AND CONCLUSIONS 

To summarise, we have confirmed that the dynamo action 
is sensitive to the nature of magnetic dissipation and demon- 
strated that magnetic reconnections (as opposed to magnetic 
diffusion) can significantly enhance the dynamo action. We 
have explored the kinematic stage of the fluctuation dynamo 
in a chaotic flow that models hydrodynamic turbulence and in 
the ABC flow, with the only magnetic dissipation mechanism 
being the reconnection of magnetic lines implemented in a di- 
rect manner. In our model, where magnetic dissipation is sup- 



pressed at all scales exceeding a certain scale do, the growth 
rate of magnetic field exceeds that of the magnetic diffusion- 
based fluctuation dynamo with the same velocity field. Even 
when the velocity field of the reconnection-based dynamo is 
reduced in magnitude as to achieve similar growth rates of 
magnetic energy density, the rate of conversion of magnetic 
energy into heat in the reconnection dynamo is an order of 
magnitude larger than in the corresponding diffusion-based 
dynamo. Thus, reconnections more efficiently convert the ki- 
netic energy of the plasma flow into heat, in our case with the 
mediation of the dynamo action. This result, here obtained 
for a kinematic dynamo, can have serious implications for the 
heating of rarefied, hot plasmas where magnetic reconnections 
dominate over magnetic diffusion (such as the corona of the 
Sun and star, galaxies and accretion discs). 

It is intriguing that reconnections play the same role lE flsll 
of converting kinetic energy into heat in superfluids and Bose- 
Einstein condensates, fluids near absolute zero at the opposite 
end of the temperature spectrum. 

Our model can be viewed as a numerical implementation 
of the elusive limiting regime of infinitely large magnetic 
Reynolds number, where magnetic dissipation can be safely 
neglected at all large scales but plays a crucial role at a cer- 
tain very small scale (we are grateful to Alex Schekochihin 
for suggesting this idea). 

In contrast to the fluctuation dynamo based on magnetic 
diffusion, the probability distribution function of the energy 
released in the flux rope dynamo has a power law form not 
dissimilar to that observed for the Solar flares. This is also 
true for the nonlinear states of the dynamo. 

The reconnection-based dynamo model suggested here can 
be generalised to include the modification of the velocity field 
by the Lorentz force. More precisely, magnetic pressure is 
assumed to be balanced by the gas pressure, so that only mag- 
netic tension needs to be explicitly included into the Navier- 
Stokes equation. Magnetic tension can readily be calculated in 
our model where magnetic field is defined only at discrete po- 
sitions of closed magnetic loops. We suggest two approxima- 
tions for the Navier-Stokes equation, one designed to model 
Alfven waves and the other suitable for the studies of nonlin- 
ear dynamos. The former model can be useful in the studies 
of nonlinear interaction of Alfven waves and Alfvenic turbu- 
lence. 

Unlike most - if not all - other simulations of the fluctua- 
tion dynamo, our computations start with a spatially localised 
initial magnetic field. This has allowed us to observe that the 
magnetised region spreads during the kinematic dynamo stage 
but its size stops growing in the nonlinear stage. This can be 
naturally interpreted as the suppression of the turbulent mag- 
netic diffusion in the saturated dynamo state. This is broadly 
equivalent to the reduction of the effective magnetic Reynolds 
number down to its marginal value (with respect to the dy- 
namo action). 

Our model of magnetic field evolution, based on tracing 
closed magnetic loops can be fruitfully applied in other nu- 
merical approaches to magnetohydrodynamics. One of well- 
known difficulties in the generalisation of smoothed-particle 
hydrodynamics to include magnetic fields is the implemen- 
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tation of the solenoidality of magnetic field. Quite notably, 
our approach satisfies the magnetic solenoidality condition 
perfectly since the modelled magnetic lines are closed at all 
times. A similar approach may be fruitful in smoothed- 
particle magnetohydrodynamics codes. 
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